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Abstract - Self-propulsion of a Janus droplet in a solution of surfactant, which reacts on a 
half of a drop surface, is studied theoretically. The droplet acts as a catalytic motor creating 
a concentration gradient, which generates its surface-tension-driven motion; the self-propulsion 
speed is rather high, 60 pm/s and more. This catalytic motor has several advantages over other 
micromotors: simple manufacturing, easily attained neutral buoyancy. In contrast to a single¬ 
fluid droplet, which demonstrates a self-propulsion as a result of symmetry breaking instability, 
for Janus one no stability threshold exists; hence, the droplet radius can be scaled down to 
micrometers. 


Introduction. — A micromotor is a particle of a mi¬ 
crometer size which is able to move progressively through 
the liquid or suspension in the absence of external fields 
and/or gradients. The motors are needed for an enor¬ 
mous number of modern applications and microtechnolo¬ 
gies, such as drug delivery, mm design of smart mate¬ 
rials, [3] etc. However, manufacturing such motors is a 
challenging task because quotidian ways of swimming are 
heavily based on inertia and, therefore, they do not op¬ 
erate at microscale (at small Reynolds numbers); alterna¬ 
tive ways, such as squirming motion [4] or ciliary propul¬ 
sion [5], widespread for living microorganisms are hardly 
practical with today’s level of technology development. 

One of the promising concepts to design a self-propelling 
particle is the so-called catalytic motor, a particle with 
nonuniform chemical reactivity along its surface, e.g. with 
a catalytic patch. Being suspended in a reactive medium, 
such a motor creates a gradient of concentration, which 
pushes the particle either diffusiophoretically or, in the 
case of charged products/reactants, electrophoretically. 
Since the particle forms the gradient of concentration 
by itself, the phenomena are also often termed as self- 
phoretic. The summary of main achievements and chal¬ 
lenges in this area can be found in the recent surveys, 
Refs. [6H9]. 

A bubble or a drop demonstrates “thermophoresis” 
(“diffusiophoresis”) being placed in a gradient of the tem¬ 
perature (concentration). The discovery of the surface- 
tension-driven migration of a drop is usually prescribed 
to Ref. [To], although it was revealed several years be¬ 


fore by Fedosov m (in Russian, for English translation 
see http://arxiv.org/abs/1303.0243vl). Dynamics of an 
interfacial drop under external temperature gradient par¬ 
allel to the interface was studied in Ref. m- 

A self-propulsion of active droplet was analyzed theo¬ 
retically in Ref. [13] and a series of other papers from the 
same group (see the detailed survey in Ref. [H]). Experi¬ 
mental observations can be found in Refs. [SHH]; unfor¬ 
tunately, no thorough comparison of the theory and ex¬ 
periments has been performed. Eor a single-fluid droplet 
the self-propulsion occurs as a result of instability. In¬ 
deed, a single-fluid droplet is spherically symmetric and 
any scalar held around it (the solute concentration, tem¬ 
perature, etc.) also obeys the spherical symmetry. In 
order to gain a self-propulsion, one has to break this sym¬ 
metry creating the radial gradient of that scalar exceeding 
a certain threshold. This makes impossible to scale down 
such a surface-tension-driven motor to the size less than 
0.1 mm. 

The natural idea is to deal with an asymmetric droplet, 
where the self-propulsion occurs without any threshold 
similarly to the so-called camphor boat [TBl - ED] (see also 
references therein). A Janus droplet (JD)—compound 
droplet with a triple line—has an asymmetry a priory and, 
therefore, it is able to move autonomously. Microfluidic 
manufacturing the JD is now a quotidian routine, [21H23j 
but usually JD is used only as an intermediate step in syn¬ 
thesis of Janus particles. Theoretically mainly the equi¬ 
librium shape of JD is studied, [24ll25] the simplest Stokes 
problem for JD is considered in Ref. [26] . 
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In this paper, we study the solutocapillary mechanism 
of self-propulsion for a perfect Janus droplet (JD), which 
comprises two attached hemispherical domains occupied 
by two different liquids. We consider a minimal model, 
which guarantees an autonomous motion, calculating the 
variation of the self-propulsion velocity with the parame¬ 
ters of liquids. Estimations of the velocity as well as the 
ways to improve the model are also discussed. 

Steady flow and diffusion in the Janus-drop sys¬ 
tem. — We consider a perfect JD, a compound drop com¬ 
posed by two hemispherical domains of a radius a occupied 
by different liquids, see Fig. [T] Generally, the surface ten¬ 
sions force the system to maintain spherical shape of all 
three interfaces; the external surfaces of the liquid domains 
1 and 2 are incomplete spheres separated by a spherically 
deformed interface. For a small internal surface tension 
coefficient compared to the external surface tension coeffi¬ 
cients, the latter two coefficients are known from chemical 
physics to be nearly equal to one another. In this case, 
the interface between liquids I and 2 is nearly flat and 
the two liquid domains are nearly hemispherical |26j . The 
case of small internal surface tension is typical enough and 
we restrict our consideration to this case. 

The drop is suspended in a solution of a surfactant, 
which is adsorbed at the external JD surface. For simplic¬ 
ity, the surfactant is assumed to be insoluble in the internal 
liquids. The simplest chemical reaction—a fixed flux of the 
reactant js —takes place at the boundary between liquids 
0 (the ambient) and 1 (reactive liquid). This assumption 
is not crucial, but it reduces the number of dimensionless 
parameters; an extension of the problem to the first order 
reaction and two species, reactant and product, can be 
carried out similarly to Ref. m- The chemical reaction 
creates a nonuniform distribution of the surfactant along 
the drop surface. 

The difference between JD and Janus particles is not 
merely hydrodynamical (see Fig. |4]for a sample flow pat¬ 
tern within the drop) but also in a continuous renewing of 
the interface between the ambient liquid and the droplet 
liquids, which diminishes the possible role of surface dif¬ 
fusion of the surfactant, allowing one to neglect it. 

Neglecting the adsorption-desorption kinetics, we as¬ 
sume that the interfacial tensions at the drop surface de¬ 
pend linearly on the surfactant concentration C at the 
external JD surface: 

7oj = 7o°^ - iCjC, j = 1 , 2 . 

In the absence of the surfactant, the interfacial tensions 
at the drop surface have to be equal to each other to deal 
with the perfect JD, [24H26] this is why is indepen¬ 
dent of j. The second terms are also small compared to 
the first one due to the same reason. Therefore, because 
of the chemical reaction at the JD surface there appears a 
gradient of the interfacial tension, which triggers a soluto¬ 
capillary motion and results in the JD self-propulsion. In 
order to simplify the analysis, we assume that the inverse 



solution of 
a surfactant (0) 


Fig. 1: (Color online) Problem geometry and coordinate sys¬ 
tem: Janus droplet comprises two hemispherical domains (liq¬ 
uids 1 and 2) suspended in a solution of a surfactant (liquid 0). 
The chemical reaction (a production of the surfactant) takes 
place at the boundary between liquids 0 and 1. 


capillary, Reynolds, and Peclet numbers associated with 
the Marangoni convection are small; this allows us to dis¬ 
regard the flow-induced interface deformations, nonlinear 
terms in the Navier-Stokes equation, and advection of the 
surfactant, respectively. The final assumption decouples 
the boundary value problem for the surfactant concentra¬ 
tion from fluid mechanics similarly to Refs. [551113] (the 
advection is taken into account in Refs. HZllSD]). In this 
limit the concentration is governed by 

V^C = 0 at r- > I, (la) 

VnC = —h(d) at r = 1 , (lb) 

C —>■ 0 at r 1 . (Ic) 


where h(d) = 1 at d < 7 r /2 and /i(d) = 0 otherwise; n 
is the external normal to the JD surface. This problem is 
non-dimensionalized using a for the lengthscale and jsO./D 
for the solute concentration {D is the surfactant diffusiv- 
ity). With (7ci-h7C2)Jsa/(2D77o) and (7C1+7C2)Js/(2D) 
as the scales for the velocity and pressure, respectively, 
the boundary value problem governing the solutocapillary 
convection around and in the JD reads: 


V • = 

0 , - = 0 , 

( 2 a) 

0 

II 

0 

[CT„r] = -ic'^rC at r = 1, 

( 2 b) 

= 0 

Is |(^^ 

II 

0 

II 

g 
, b 

II 

( 2 c) 


= -UspBz at r » 1 . 

(2d) 


Here 7 c = 7 ci at d < 7 r /2 and 7 c = 702 otherwise; the 
brackets denote the jump of the corresponding quantity 
across the drop surface: [/] = /i ^2 — /o- The sought veloc¬ 
ity of self-propulsion Usp is determined by the condition 
of no overall force imposed to JD. In other words, this 
velocity is chosen to compensate the Marangoni-induced 
drag. 

Four dimensionless parameters enter the problem 
boundary value problem; they are 


27cj 

I 3 = 7Cj = ^ 

7o 7C1 + 7C2 


(J = l,2), 
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Fig. 2: (Color online) The distribution of surfactant concen¬ 
tration according to bright region corresponds to higher 
C, dark one is for the lower C. 


which represent the dimensionless measures of the inter¬ 
nal viscosities and solutocapillary effect, respectively. By 
definition, the second pair of parameters is coupled by the 
relation yci + 7 C 2 = 2, therefore, Ugp is governed by three 
independent dimensionless parameters only. 

The solution for the concentration is shown in Fig. [2l it 
is given by [28l[29] : 


1 CnP2n+litJ‘) 

2 j - ' ' J ^ 2 n -|-2 

n =0 


Cr, 


4n -b 3 
8 (n-b 1)3 


^"2n(0), (3) 


where fj, = cos'd and Pn are the Legendre polynomials. 
m The concentration is higher near the reactive surface, 
decreasing as -d grows beyond 7 r/ 2 . At large distance from 
the JD, the monopole contribution l/(2r) prevails. 

The solution for the velocities can be presented in terms 
of the vector potentials i^(l> is th® unit vector for the 
azimuthal angle (j)) introduced in such a way that Vj = 
V X Therefore, the amplitude of vector potential 

Ip is related to the conventional streamfunction ipc |32j via 
tpc = rsinOip. 

For the ambient ip can be presented as 


n—1 


-^n^nl (/^) 


(1 - r^) -b Usptpoo, 




1 

4 



sind, 


(4) 

(5) 


where P„i(^) = ^/l^^y?dPn{^j)/d^i are the associate Leg¬ 
endre polynomials. m Since Ugp is the velocity of self¬ 
propulsion (to be determined), the term ipoo represents 
the Stokes flow past a solid particle of unit radius; any 
other field, which tends to — ^rsin'd far from the drop is 
also appropriate. For the internal fluids the expressions 
found in Ref. |26j can be implemented. 

The ansatz for the velocity fields satisfies the boundary 
conditions at the internal interface; demanding the rest of 
boundary conditions, we end up with the linear algebraic 
set for the coefficients and appropriate coefficients for 
This set of equations is solved numerically with the 
Maple software applied. 



[ri]/(2ri) 


Fig. 3: (Color online) The variation of the velocity of self¬ 
propulsion Uap (in units [pci + PC 2 )jsa/ 2 Dfjo) with the 
weighted viscosity difference [ri]/{2ri) = (572 — + 72 ) for 

three values of the mean viscosity of the JD r] = (r)i-bj)2)/(277o): 
7 = 0.1 (lines 1), 7 = 1 (lines 2), and 77 = 10 (lines 3); the solid 
lines correspond to [ 7 c] = ( 7 C 2 — 7 ci)/(' 7 ci + ■ 702 ) = 0.5; the 
dashed lines are for [ 7 c] = 0. In this case Uap is an even func¬ 
tion of 772 — 71 • In the inset the zoomed in lines 1 are shown. 



Fig. 4: Streamlines (isolines of the streamfunction ipc = 
rsin'di/j) for rji — 0.5, 72 = 1.5; panel (a): 701 = 7 C 2 = 1; 
panel (b): 701 = 1.5, 702 = 0.5. The difference between the 
streamlines is 0.002 for the outer flow and 0.0002 for the flow 
in the JD. 


The force imposed to the JD is given by [33] 2tt{3Usp — 
4Ai); setting this value to zero one can determine the ve¬ 
locity of self-propulsion Ugp- [Note that Ai also contains 
the contribution linear in Ugp — the overall terms propor¬ 
tional to Ugp is the resistance (inverse mobility) found in 
Ref. [33].] 

Results and Discussion. — We start the discussion 
with the case 7 ci = 702 = 1 , see the dashed lines in 
Fig. [31 In this case the velocity of self-propulsion is deter¬ 
mined by 71.2 only or by their combinations 77 = (r]i+r] 2)/2 
and [ 7 ] = 72 — 7 i, which are the mean JD viscosity and 
the difference of the internal viscosities, respectively. The 
velocity of self-propulsion is even function of [ 7 ] and it 
reaches the maximum at larger contrast of the viscosities, 
72 7 i- It means that the intensive motion in the first 

(less viscous) internal liquid guarantees the higher values 
of Ugp, even though the convection in the second internal 
liquid is damped. This trend, however, is violated at small 
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rj, see the inset in Fig. [31 when the JD is almost inviscid 
and the distribution of the small mean viscosity between 
the internal liquids becomes of a little importance. In 
other words, the motion in each liquid and Ugp are almost 
independent of [rj\. The self-propulsion velocity is affected 
by subtle details and Ugp becomes nonmonotonic func¬ 
tion of [r;]: it initially decreases at small [ry], approaches a 
minimum, and then grows in compliance with the general 
tendency. 

If the solutocapillary effects at the external interfaces 
are different, the symmetry of Uap with respect to inver¬ 
sion of [ 77 ] is lost, see Fig. [31 but Usp is still invariant 
under the transformation of 771 -H- 772 and yci 7 C 2 (at 
fixed direction of self-propulsion). As one can intuitively 
expect, the larger velocity of self-propulsion is attained if 
the liquid with larger 7 c (more intensive generation of the 
fluid motion) has smaller viscosity (smaller hindrance to 
the motion). Again, for smaller 77 , the difference in the 
internal viscosities is less important. 

The flow structure is shown in Fig. |4l Qualitatively the 
flow looks similar to that for the Stokes flow past JD, [2^ 
a toroidal vortex in each internal fluid. It should be em¬ 
phasized that, similarly to the cited paper, at the external 
boundary of the JD the internal fluids are comoving, but 
they move oppositely near the internal interface. This 
leads to formation of a small-intensity vortex with the op¬ 
posite direction of fluid rotation adjacent to the internal 
interface, which provides matching of the velocity there. 
[This vortex is usually situated in a liquid with smaller 
intensity of fluid motion, but the opposite is also possible 
under certain conditions, see Fig. SJa).] The typical ve¬ 
locity is by several orders of magnitude smaller there than 
for the rest of fluid, therefore, on the scale of the figure 
one can see only the boundary of this vortex (the line cor¬ 
responding to tp = 0), the streamline which starts at the 
internal interface and ends up at the z-axis. This region 
hence can be thought of as a stagnant zone. 

Conclusions. — A self-propulsion of a JD in a so¬ 
lution of surfactant, which experiences the zeroth order 
chemical reaction (fixed flux) at one of the external in¬ 
terfaces of JD. Using the minimal model, we derive the 
dimensionless velocity of self-propulsion as a function of 
the internal viscosities and the solutocapillary constants. 
It is worth noting that most of the restrictions (a per¬ 
fect JD, the fast sorption kinetics, oversimplified chemical 
reaction, insolubility of the surfactant in both internal liq¬ 
uids, etc.) can be relaxed, although the corresponding 
generalizations would mainly need a numerical analysis. 

In order to estimate the typical velocity, we consider 
the JD of I /im radius in a castor oil 770 « 10 P 
(IPa-s); the typical variation of the interfacial tension 
( 7 C 1 + 7C2)jsa/(2D) is 1 dyne/cm, which corresponds to 
a 1% variation of the surface tension along the surface. (In 
fact, the estimation below does not depend on the radius 
at given difference of the surface tension.) For Usp = 0.06 
(?7 can be easily made small for so viscous ambient), the di¬ 


mensional velocity of self-propulsion is as large as 60 /rm/s. 
This value is rather large and it can be further augmented 
by choosing the less viscous ambient. 

Among the main advantages of this kind of catalytic 
motor, the following can be pointed out: (i) simple man¬ 
ufacturing via the well-developed techniques; (ii) simple 
suspending JD in a liquid, just varying the densities of 
internal fluids (by now most experiments are performed 
near the solid substrate with the only exception. Ref. [33]); 
(iii) sufficiently high velocity of self-propulsion. 

* * * 
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